Optical conductivity of a superfluid density wave 
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We present a calculation of the low frequency optical conductivity of a superconductor in the 
presence of quenched inhomogeneity in both the superfluid and normal fluid densities. We find that 
inhomogeneity in the superfluid density displaces spectral weight from the condensate to a frequency 
range that depends critically on the spatial correlation of normal and superfluid density fluctuations. 
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Measurements of the optical conductivity, cr{uj), in the 
microwave and terahertz regimes have proved to be a 
valuable probe of quasiparticle dynamics in cuprate su- 
perconductors. Because even relatively poor samples are 
in the clean limit where the scattering rate is smaller 
than the gap frequency, or 1/r ^ A, the spectral weight 
of the dissipation due to quasiparticle scattering is far 
greater than that associated with pair creation. As a 
consequence the quasiparticle Drude peak can be clearly 
resolved, and its width can yield the scattering rate. Im- 
portant clues as to the nature of the scattering can be 
obtained by comparing the rates measured by microwave 
conductivity with those measured by angle-resolved pho- 
toemission (ARPES). 

Extracting 1/t from the microwave conductivity of a 
cuprate requires ct^lu) to be consistent with the two-fluid 
model for a clean, d-wave superconductor In this 
model the total conductivity contains two contributions: 
a (5-function ed lu — due to the condensate and a Drude- 
like peak due to quasiparticle scattering. The spectral 
weights of these components are the super and normal 
fluid densities, respectively. As the temperature T tends 
to zero the normal fluid density, p„, is expected to de- 
crease linearly with T. The superfluid density, ps is ex- 
pected to increase at the same rate, so that the total 
density remains constant. This simple behavior permits 
Pn.s and 1/t to be determined in the superconducting 
state. 

The two-fluid model provides an excellent description 
of the microwave properties of Y Ba2Cu2,0'j-s (YBCO) 
single crystals . The scattering rates inferred from the 
two-fluid analysis are remarkably small, approaching ~ 
10 GHz as T -> 0. On the other hand, analysis of ARPES 
yields low temperature scattering rates that are larger by 
a factor of at least 100 |^. One possible explanation is 
that the single particle self-energy measured by ARPES 
is very different from the momentum relaxation rate that 
determines the conductivity. However, an alternative ex- 
planation is that the ARPES data are collected almost 
exclusively from studies of Bi2Sr2CaCu20s+s (BSCCO) 
rather than YBCO. 

It would seem that the microwave data on BSCCO 



would hold the resolution to this question. However, 
analysis of such data has been problematic because the 
spectra are not consistent with the two-fluid model de- 
scribed previously. In fact the success of the two-fluid 
description is unique to YBCO single crystals - all other 
cuprates that have been measured show distinctly dif- 
ferent properties In these materials the spectral 
weight of the Drude peak extrapolates to a large residual 
value as T ^ 0, suggesting that a substantial fraction of 
the normal fluid fails to condense. 

Recently, the frequency dependence of the conductivity 
was reported in optimally doped BSCCO [||. The mea- 
surements were performed using coherent time-domain 
spectroscopy, a technique which can determine cr(w) up 
to frequencies uj/27t ~1 THz. The width of the Drude 
peak at low temperature was determined to be '^300 
GHz. Integration of the real part of the conductiv- 
ity over the experimental range showed that the uncon- 
densed portion is at least 30% of the condensate spectral 
weight. A study of BSCCO for several carrier concentra- 
tions showed that this percentage increases with doping, 
reaching about 60% in an overdoped sample 

The two-fluid description of (j(uj, T) in such samples is 
not internally consistent. Fitting the spectra requires p„ 
to vary with frequency and actually increase with de- 
creasing T for ui <300GHz. On the other hand, the 
quasiparticle conductivity could be described by phys- 
ically reasonable parameters if Re a has a third com- 
ponent: a peak centered at uj—Q, whose spectral weight 
grows in proportion to ps as T decreases. 

The proportionality to ps suggests that this anomalous 
part of a is related to current fluctuations of the super- 
fluid. However, in a homogeneous superconductor the 
spectral weight of such fluctuations decreases exponen- 
tially with T, which is opposite to the growth of spectral 
weight that was observed. Recently, Barabash et al. |Tc| ] 
showed that a component of Re a can increase as T — > if 
Ps is spatially inhomogeneous. They considered a Joseph- 
son junction network with a quenched variation 5 J about 
the mean coupling J. (J is the superfluid density defined 
on a lattice rather than a continuum). They showed that 
the global phase stiffness J equals J — {5,P)/ J. J is less 
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than J because the phase varies more rapidly in regions 
where the stiffness is below the average value. 

The result stated above is important to the optical con- 
ductivity because the spectral weight of the condensate 
i5-function is proportional to J, whereas the total spectral 
weight of condensate current fluctuations is proportional 
to J. Thus inhomogeneity of the superfluid density must 
displace spectral weight ^ {5,P)/J from the (5-function 
to nonzero frequency, li 5J varies in proportion to J 
then the displaced spectral weight will track that of the 
condensate, in agreement with terahertz experiments Q . 

The existence of strong spatial inhomogeneity in 
BSCCO has been demonstrated by scanning tunneling 
microscopy (STM) measurements The local den- 

sity of states (LDOS) of BSCCO varies randomly in 
space, with spatial fluctuations that have a minimum 
wavelength of '^50 A. The variations in LDOS suggest 
quenched inhomogeneity in local carrier concentration, 
X, and therefore in the local ps- 

To determine whether quenched inhomogeneity in ps 
is responsible for the anomaly in cr(Lj), we need to con- 
sider the spectrum of the spectral weight removed from 
the condensate. In a system with no normal component 
the spectral weight would be expected to appear at the 
natural oscillation frequency of the condensate, which is 
the Josephson plasma frequency, Ug ■ In optimal cuprates 
cjs/27r^200 THz, whereas the anomalous dissipation is 
found below '^1 THz. Inhomogeneity can only explain 
the data if the presence of a normal fluid in addition 
to the superfluid causes the displaced spectral weight to 
appear at frequencies much smaller than ujs- 

In this paper we calculate the change in a{uj) due to 
inhomogeneity in a superconductor in which normal and 
superfluid coexist. This is the flrst calculation in which 
the spectral weight of the conductivity is consistent with 
the theorem of Ref. [0. We flnd that (j{ui) is extremely 
sensitive to correlations of the normal and superfluid den- 
sity variations in the the medium. The spectrum depends 
crucially on whether the correlation is positive (such that 
regions of large ps have large p„) or negative. If the cor- 
relation is positive, the displaced spectral weight indeed 
shifts to the plasma frequency. However, when the fluc- 
tuations of normal and superfluid density are anticorre- 
lated, the spectral weight appears at low frequencies, in 
agreement with microwave and terahertz measurements 
on cuprate superconductors. 

To treat the conductivity in the presence of inhomo- 
geneity we apply the extended two-fluid phenomenol- 
ogy developed by Pethick and Smith jl^], and Kadin 
and Goldman [|l3|. This approach successfully describes 
quenched inhomogeneity at the normal-superconductor 
interface and fluctuating inhomogeneity, as in the 
Carlson- Goldman oscillations iQ. In the extended two- 
fluid model the superfluid is accelerated by gradients of 
the chemical potential as well as electric fields, that is. 



Js = Ps{E-V^J^s)■ (1) 

The chemical potential has the subscript s because in 
a superconductor ji is the energy per electron required 
to add a pair to the condensate. The corresponding 
equation for the normal fluid current requires solving 
the Boltzmann equation for the quasiparticle distribution 
function. However for frequencies less than 1/t, the dis- 
tribution function is the equilibrium distribution shifted 
by the 'quasiparticle chemical potential' or If the 
normal fluid is in local equilibrium with the condensate 
Pn — ^J■s, which differs from "global" equilibrium where 
/i„=0. In the low frequency regime the constitutive rela- 
tion for the normal fluid has the simple form: 

Jn^ Pn{E -Wpn) - Jn/r. (2) 

A closed system of equations requires continuity rela- 
tions. The total charge of the superconductor separates 
naturally into a normal component, Qn that depends on 
both the coherence factors and the distribution function, 

<9« = ^qkfk (3) 

k 

where = uj. — vj,, and a superfluid component that 
depends only on coherence factors, 

Qs = ^2evl = 2eNFPs. (4) 

k 

Under conditions for which /i„ can be defined, the normal 
fluid charge is given by, 

Qn = 2Nf\{Ps - Mn), (5) 

where Np is the normal state density of states at the 
Fermi level. The parameter A relates the normal fiuid 
charge to the shift of pn away from local equilibrium. 

In a superconducting medium the normal and super- 
fluid charge are not separately conserved. Intercon- 
version of normal and superconducting charge occurs 
through two types of processes. In the flrst process, Qn 
changes as quasiparticles recombine or scatter. In the 
second, the quasiparticle charge changes even if fk re- 
mains constant. In this process, Qn varies because the 
quasiparticle excitation spectrum, and consequently the 
effective charge, adjusts to the local value of ps- Conti- 
nuity equations that include both types of exchange be- 
tween the two fluids are: 

Qn.s+V ■Jn,s = {-,+){— -\Qs), (6) 

where tq is the rate of conversion of normal charge into 
superfluid charge due to scattering and recombination 
processes. The above system of equations is closed by 
V-E = Q/eo. 
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To see how quenched inhomogeneity affects the optical 
conductivity, we consider the simplest possible model: a 
one-dimensional sinusoidal variation in normal and su- 
perfluid density. We assume that the densities of the two 
components vary as Ps,n{x) = {ps,n) + ^e{psq,„qe*«^}, 
where (p) is the average density and pq is the ampli- 
tude of the inhomogeneity at wavevector q. Because the 
medium is inhomogeneous, a uniform applied field gener- 
ates a field at q. Solving the extended two- fluid equations 
to lowest order in Pnq,sq, we obtain, 



Eo 



Psq Pnq^ 



where, 



.2-c.2(l-A) + K-zc.A/r)f^' 



^{uj + i/T*Q)-vlq^l{l-\) 
{LO + i/T){u + i/T*Q)~vlq^/y 



(7) 



(8) 



= Ps^n/^Q + w^.„q^ = Ps^n / '2.N pe^ , and r* = 
tq{1-\). 

The uniform current density in response to these fields 
is given by Jq = cto-E'o + <^qE-q, where uq is the uniform 
two-fluid conductivity, and Gq is the conductivity that 
varies with wavevector q. In the two fluid model these 
are given by, 



^Ps^sq . ^Pn. 



nq 



(9) 



The effective conductivity of the medium, a = ai + ia2, 
is the ratio of the uniform current density to the uniform 
field, so that, a = aa + aqE-q/ Eq. The second term in 
this equation is the change in the optical conductivity 
due to inhomogeneity, or Act. 

The extra term in the conductivity is particularly sim- 
ple if p„ = 0, in which case. 



Act2 = - 



Psq Psq 



(10) 



The inhomogeneity in the superfiuid density indeed re- 
moves spectral weight p\qlps from the condensate 5- 
function, in agreement with the results of Ref . . In the 
absence of a normal fluid component the spectral weight 
reappears in a J-function at the Josephson plasma fre- 
quency. 

We next assume that Pn 7^ (even at T = 0), and 
describe how this assumption affects Act(w). We focus 
on the behavior of Act when the normal fluid density 
fluctuations are either perfectly correlated or anticorre- 
lated with those of the superfluid density. We take for 
two-fluid parameters values that are suggested by the 
terahertz and microwave experiments: (ps/eo)"'^^^=1000 
THz, (p„/eo)^/^=800 THz, and r"! = Tq1=3 THz. The 
existence of normal fluid at T = implies a large den- 
sity of states at the Fermi level, Nq. While the origin of 
iVo is outside the scope of this paper, one may speculate 



that it is closely related to the strong inhomogeneity ob- 
served by STM. If we make the reasonable approximation 
of neglecting BCS coherence factors for the quasiparticle 
states introduced by disorder, then A = No/Np. 
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FIG. 1. AfTi as a function of frequency (tj/27r) for anticor- 
related variations in ps and p„. Spectral weight decreases for 
increasing Vaqr: 0.25, 0.5, 1.0, and 2.0. The same curves are 
shown in a normalized plot in the inset. 

We begin with the case where the density fluctuations 
are perfectly anticorrelated, so that p^q = —pnq and the 
total fluid density is uniform throughout the medium. 
Fig. 1 shows Acti(w) with {psq/eo)^^'^=100 THz, for sev- 
eral values of Vsqr. Acti(cj) is positive and centered at 
Lu = rather than lUs- The spectra depend strongly on 
Vsqr. For Vgqr <d, Acti has a Drude-like spectrum, 
whose width is ~ 1/r. As Vsqr increases beyond unity, 
the spectral weight drops and a peak near the Carlson- 
Goldman frequency Vsq appears in the spectrum. 

The key issue is the fraction of the displaced conden- 
sate spectral weight that appears in the low-frequency 
peak, as opposed to w ~ w^. In Fig. 2 we compare the re- 
duction in condensate weight with the increase in dissipa- 
tion at low frequency. The change in condensate spectral 
weight was evaluated from the limi^^o{TT/2)LL!Aai{Lu). 
The low frequency spectral weight was obtained by nu- 
merically computing the integral of Acti with respect to 
uj from to 100 THz. Fig. 2 shows these two quanti- 
ties, normalized to p'^^/ps, as a function of Vsqr. They 
are equal in magnitude but opposite in sign, which shows 
that all of the spectral weight removed from the conden- 
sate appears at low frequency and none appears at cus- 
Moreover, the decrease of condensate spectral weight co- 
incides exactly with the prediction of Barabash et al. |ic| ] 
as Vsqr — > 0, but vanishes for Vgqr 3>1. 

The surprising results presented above are a straight- 
forward consequence of the ant icor relation of the density 
fluctuations. There is no dissipation near oj^ or ujp = 
[{ps -f Pri)/eo]^^^ because the inhomogeneous medium 
is dynamically homogeneous at high frequencies. For 
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FIG. 2. Comparison of the spectral weight displaced from 
the condensate by the quenched variation in pa and the spec- 
tral weight that appears in Acri at low frequency. 



a; ~ the fractional difference in the conductivity of 
the superfluid rich and superfluid poor regions is only 
~ l/(wpT)^. The anomalous dissipation appears instead 
at low frequency where the conductivity is strongly inho- 
mogeneous. a is smaller in the regions that are superfluid 
poor and normal fluid rich. E will be larger in such re- 
gions, which is precisely equivalent to more rapid phase 
variation in regions of low stiffness. Thus the increase 
in dissipation arises ultimately from an amplification of 
the electric field in regions with greater than average p„. 
Finally, the dynamical inhomogoncity disappears when 
ij) ^ Vs(l because the normal and superfluid response 
again become indistinguishable in this regime. 
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FIG. 3. Aai vs. frequency for the same parameters as Fig. 
1, except variations in ps and pn are positively correlated. 

We are now prepared to understand the sharp differ- 
ence in Act if the super and normal fluid density waves are 
in phase. Fig. 3 shows Re{A(7{uj)} calculated with iden- 
tical parameters as in Fig. 1, except that Ps = Pn- The 
change in conductivity is negative, and the spectra are 



nearly independent of q in contrast to the strongly q de- 
pendent and positive change generated when ps = —pn- 
The reduction in conductivity is exactly as expected from 
the previous arguments: now regions of large E coin- 
cide with normal fluid poor regions and the dissipation 
is attenuated. The spectra are nearly independent of q 
because the dynamical inhomogoncity docs not tend to 
zero when the normal and superfluid response become in- 
distinguishable. The spectral weight is displaced equally 
from the condensate and the normal fluid, and weight 
2plg/ps shifts to LOp. 

We conclude with a qualitative comparison of the cal- 
culation and the low frequency a{uj) of the cupratcs. De- 
spite an oversimplified representation of quenched dis- 
order, the model provides a natural explanation for the 
most surprising feature of the data, the dramatic increase 
of the anomalous spectral weight as x increases from un- 
der to overdoped values. The calculation shows that this 
can occur if the spatial correlation of ps and /?„ depends 
strongly on x. Such a change in the nature of the corre- 
lations upon crossing optimal doping Xopt would indeed 
be expected, based on very general considerations. For 
a sample with x < Xopt, regions of lower than average 
Ps will be strongly underdoped patches. These patches 
will be insulators, or at least regions with small On 
the other hand, the regions of low ps in a sample with 
X > Xopt will be strongly overdoped, i.e. regions of large 
Pn- Thus the correlation varies systematically from posi- 
tive to negative with increasing x, leading to a dramatic 
increase in the low frequency spectral weight. 
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